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Two-band Fluctuation Exchange Study on the Superconductivity of 
/?'- (BEDT-TTF) 2 IC1 2 under High Pressure. 

Tsuguhito Nakano * and Kazuhiko Kuroki 

Department of Applied Physics and Chemistry, The University of Electro- Communications, Chofu, Tokyo 182-8585, Japan 

We study the pressure dependence of the superconducting transition temperature of an or- 
ganic superconductor /3'-(BEDT-TTF)2lCl2 by applying the fluctuation exchange method to 
the Hubbard model on the original two-band lattice at 3/4- filling rather than the single band 
model in the strong dimerization limit. Our study is motivated by the fact that hopping param- 
eters evaluated from a first-principles study suggest that the dimerization of the BEDT-TTF 
molecules is not so strong especially at high pressure. Solving the linearized Eliashberg's equa- 
tion, a d X!/ -wave-like superconducting state with realistic values of T c is obtained in a pressure 
regime somewhat higher than the actual experimental result. These results are similar to those 
obtained within the single band model in the previous study by Kino et al. We conclude that 
the resemblance to the dimer limit is due to a combination of a good Fermi surface nesting, a 
large density of states near the Fermi level, and a moderate dimerization, which cooperatively 
enhance electron correlation effects and also the superconducting T c . 

KEYWORDS: superconductivity,/3'-(BEDT-TTF) 2 ICI 2 , pressure effects.multi band, FLEX 
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1. Introduction 

Over the past few decades, a considerable amount of 
studies have been performed on organic conductors both 
experimentally and theoretically. Especially, the occur- 
rence of unconventional superconductivity in these ma- 
terials, where low dimensionality and/or strong electron 
correlation may be playing an important role, has be- 
come one of the fascinating issues in condensed matter 
physics. 1 

The title material of this paper, /?'-(BEDT-TTF) 2 ICl 2 
is a charge transfer organic material which consists of 
cation BEDT-TTF (abbreviated as ET) molecule layers 
and anion ICI2 layers. This material is a paramagnetic 
insulator at room temperature and ambient pressure, and 
becomes an antiferromagnetic (AF) insulator below the 
Neel temperature Tjv = 22 K. Regarding the electronic 
structure, since two ET molecules are packed in a unit 
cell as shown in Fig. 1(a) with 0.5 holes per ET molecule, 
it is a 3 / 4-filled two-band system. Moreover, in the j3'- 
type arrangement of the ET molecules, which is rather 
modified from the /3-type because of the small size of the 
anion ICI2, ET molecules form dimers in the pi direction, 
which opens up a gap between the two bands. Thus, only 
the anti-bonding band intersects the Fermi level, so that 
it may be possible to look at the the system as a half-filled 
single band system. At ambient pressure, this picture is 
supported by the fact that system becomes an insulator 
despite the band being 3/4-filled. 

Recently, superconductivity has been found in this ma- 
terial under high pressure (above 8.2 GPa) by Taniguchi 
et al. It has the highest transition temperature T c 
(=T c roax = 14.2 K at p = 8.2 GPa ) among all the 
molecular charge-transfer salts. 2 Since the superconduct- 
ing phase seems to sit next to the antiferromagnetic insu- 
lating phase in the pressure-temperature phase diagram, 
there is a possibility that the pairing is due to AF spin 




Fig. 1. Schematic illustration of the ET molecule layer, (a) The 
original two band lattice. A small oval represents an ET molecule. 
pi, p2, ■ ■ ■ stand for the hopping integrals t(pl), t(p2), ■ ■ ■ . The 
shaded portion denotes the unit cell, (b) The effective single band 
lattice in the dimer limit with effective hoppings to, t\, and t 2 . 



fluctuations. In fact, Kino et al. have calculated T c and 
the Neel temperature Tjy using the fluctuation exchange 
(FLEX) 3 method on an effective single-band Hubbard 
model at 1/2-filling, namely the 'dimer model' obtained 
in the strong dimerization limit (Fig. 1(b)). 4,5 In their 
study, the hopping parameters of the original two band 
lattice are determined by fitting the tight binding disper- 
sion to those obtained from first principles calculation, 6 
and the hopping parameters of the effective one band 
lattice are obtained from a certain transformation. The 
obtained phase diagram is qualitatively similar to the ex- 
perimental one although the superconducting phase ap- 
pears in a higher pressure regime. 

Nevertheless, we believe that it is necessary to revisit 
this problem using the original two-band lattice due to 
the following reasons. (i)If we look into the values of the 
hopping integrals of the original two-band lattice, the 
dimerization is not so strong, and in fact the gap between 
the bonding and the antibonding bands is only 10 percent 
of the total band width at 8 GPa. (ii) The effective on- 
site repulsion in the dimer model is a function of hopping 
integral and thus should also be a function of pressure. 7 
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(iii) It has been known that assuming the dimerization 
limit can result in crucial problems as seen in the studies 
of re-(ET) 2 X, in which the dimer model gives d x 2_ y 2 -wave 
pairing with a moderate T c , while in the original four- 
band lattice, d x 2_ y 2 and d xy -wa.ve pairings are nearly 
degenerate with, if any, a very low T c . 8 ' 9 Note that in 
the case of k-(BEDT-TTF) 2 X, the band gap between 
the bonding and the antibonding band is more than 20 
% of the total band width, 10 which is larger than that in 
the title compound. 

In the present paper, we calculate T c and the gap func- 
tion by applying the two-band version of FLEX for the 
Hubbard model on the two-band lattice at 3/4-nlling us- 
ing the hopping parameters determined by Miyazaki et 
al. We obtain finite values of T c in a pressure regime sim- 
ilar to those in the single band approach. The present 
situation is in sharp contrast with the case of k-(ET) 2 X 
in that moderate values of T c (or we should say "high 
T c " in the sense mentioned in §4.4) are obtained even in 
a 3/4- filled system, where electron correlation effects are, 
naively speaking, not expected to be strong compared to 
true half- filled systems. The present study suggests that 
the coexistence of a good Fermi surface nesting, a large 
density of states and a moderate (not so weak) dimeriza- 
tion cooperatively enhances electron correlation effects 
and leads to results similar to those in the dimer limit. 
We conclude that these factors that enhance correlation 
effects should also be the very origin of the high T c itself 
of the title material. 

2. Formulation 

In the present study, we adopt a standard Hubbard 
Hamiltonian having two sites in a unit cell, where each 
site corresponds to an ET molecule. The kinetic energy 
part of the Hamiltonian is written as 

Wkin - £ *( C )( C (i.,i,+l),* C (W,).* 



+ 



U (*x,ij, + l),o- 



+ *(« 2 )4..i»)^ C ( < - i «). ff 

+ t(j>2)d\ ix _ liy)ta c {ix , iy)ta 
+ h.c — fi 



(1) 



where Cj j(T and di yCr are annihilation operators of electrons 
with spin a at the two different sites in the i-th unit 
cell, and \i represents chemical potential. i(pl), t(p2), ■ ■ ■ 
are the hopping parameters in the pi, p2, ■ ■ ■ directions, 
respectively. 

The interaction part is 



(2) 



where U is the on-site electron-electron interaction, and 
n i.a = A,a c i,<y and nf a = dj |0 .dj )<7 arc the number oper- 



ators. The pressure effect on the electronic structure is 
introduced through the hopping parameters within this 
model. As we have mentioned above, we use the hop- 
ping parameters determined by Miyazaki et al. 6 as shown 
in Tabled, which well reproduce the results of the first 
principles calculation. For pressures higher than 12 GPa, 
which is the highest pressure where the first principles 
calculation have been carried out, the values of the hop- 
ping integrals are obtained by linear extrapolation as in 
the previous study. 

In the present study, we have employed the two-band 
version of the fluctuation exchange (FLEX) approxima- 
tion to obtain the Green's function and the normal self- 
energy. For later discussions, let us briefly review the 
FLEX method, which is a kind of self-consistent random 
phase approximation (RPA) . Since FLEX can take large 
spin fluctuations into account, these methods have been 
applied to the studies of high-T c cuprates and other or- 
ganic superconductors. 

The (renormalized) thermal Green's function G(k,e) 
is given by the Dyson's equation, 



G~\k, e n ) = G \k, e n ) - S(fe, e n ), 



(3) 



where e n = (2n + l)nT is the Matsubara frequency with 
n = 0, ±1, ±2, • • • . Go(k, e n ) is the unperturbed thermal 
Green's function and S(fe, e n ) is the normal self-energy, 
which has an effect of suppressing T c . 

Using G(k,s n ) obtained by solving eq.(3), the irre- 
ducible susceptibility xo(<7,w m ) is given as 



Xo(q, w m ) = -jj £ G(k + q,iv m + e n )G(k, e n ) 



(4) 



k,n 



where oj m is the Matsubara frequency for bosons with 
m = 0, ±1,±2, ••• and N is the number of A;-point 
meshes. By collecting RPA-type diagrams, the effective 
interaction and the singlet pairing interaction 
are obtained as 



3 1 

= —^u 2 Xs(q,^m) - ^u 2 Xc{q,u m ) (5) 

VW(q,u m ) = U +^U 2 Xs(q^ m )-^U 2 X c(q^J$) 

where Xs,Xc are spin and charge susceptibilities, respec- 
tively, given as 

XsM^m ) = t— 77 — j r- 7 

iTUxo(q,u m ) 

Then the normal self-energy is given by 

S(fc,e„) = -^£G(fe- 9 , £ „) [vW{q,w m ) - U 2 X o(q^ m ) 

q.m 

(8) 

The obtained self-energy E(fe, e n ) is fed back into the 
Dyson's equation eq.(3), and by repeating these proce- 
dures, the self-consistent G(k,e n ) is obtained. 

In the two-band version of FLEX, G{k 1 e n ), Xo, Xs,c, 
E(fc,w m ) become 2x2 matrices, e.g. G a p, where a and 
(3 denote one of the two sites in a unit cell. 

(2) 

Once G a p(k,e n ) and are obtained by FLEX, we 
can calculate T c by solving the linearized Eliashberg's 
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Table I. Pressure dependence of the hopping parameters of 
/?'-(ET) 2 ICl 2 determined by Miyazaki and Kino, (ref.6) 
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Fig. 2. Transition temperature T c as functions of pressure for sev- 
eral values of U. 
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Fig. 3. Contour plots of |G(fc)| 2 and the nodes of the supercon- 
ducting gap function </>(k) for U = 0.9 eV at p = 15.25, 15.5, 
16.0, and 16.2 GPa. At p = 15.25 GPa, 4>(k)4>(k + Q) is positive 
at the circled portions of the Fermi surface. 



equation as follows, 

rp 

\<f> a p(k,£ n ) = -— V^{k-k',e n -e m ) 

k' ,m,a' ,/3' 

x G aa > (k', e m )Gpp> {—k', -e m )(t> a >p> (fc', e m ), (9) 

where <j>(k) is the superconducting gap function. The 
transition temperature T c is determined as the tempera- 
ture where the eigenvalue A reaches unity. 

In the actual calculation, we use 64 x 64 fc-point 
meshes and 16384 Matsubara frequencies in order to 
ensure convergence at the lowest temperature studied 
(T/\t(pl)\ = 0.002). The bandfhling (the number of elec- 
trons per site) is fixed at n = 1.5. 

When U~xo(q,Um = 0) = 1, the spin susceptibility 
diverges and a magnetic ordering takes place. In the 
FLEX calculation in two-dimensional systems, Mermin- 
Wagner's theorem is satisfied, 11 ' 12 so that Uxo(lt^m = 
0) < 1, namely true magnetic ordering does not take 
place. However, this is an artifact of adopting a purely 
two-dimensional model, while the actual material is quasi 
two dimensional. Thus, in the present study, we assume 
that if there were a weak three dimensionality, a mag- 
netic ordering with wave vector q would occur when 

max {Uxo (<?, u n = 0) > 0.995}, (10) 

i 

is satisfied in the temperature range where A < 1. There- 
fore we do not calculate T c in such a parameter regime. 13 

3. Results 

Now we move on to the results. Figure 2 shows the 
pressure dependence of T c obtained for several values 
of U. Since our calculation is restricted to tempera- 
tures above ~ 6 K, T c is obtained within that tem- 
perature range. At pressure lower than the supercon- 



ducting regime, the system is in the AF phase in the 
sense we mentioned in §2. The maximum T c obtained is 
T c max = 8.7 K (at 15.5 GPa for U = 0.9 eV, and 16.0 
GPa for U = 1.0 eV), which is somewhat smaller than 
the experimental maximum value of T c , but can be con- 
sidered as fairly realistic. The overall phase diagram is 
qualitatively similar to the experimental phase diagram, 
but the pressure range in which superconductivity oc- 
curs is above ~ 14 GPa and extends up to ~ 17 GPa 
or higher, is higher than the experimental results. These 
results are similar to those obtained within the dimcr 
model approach. 4 

Figure 3 shows the nodal lines of <fi(k) and the con- 
tour plots of \G(k)\ 2 for several values of pressure with 
U = 0.9 eV. In the plots of |G(fc)| 2 , the center of the 
densely bundled contour lines correspond to the ridges 
of |G(fc)| 2 and thus the Fermi surface, while the thick- 
ness of these bundles can be considered as a measure 
for the density of states near the Fermi level, namely, 
the thicker these bundles, the larger number of states lie 
near the Fermi level. With increasing pressure, the Fermi 
surface changes its topology from a one dimensional one 
open in the k c direction to a closed two dimensional one 
around (0, 7r). Again like in the single band approach, 
the pairing symmetry is d^y-wave-like in the sense that 
4>(k) changes its sign as (H 1 — ) along the Fermi sur- 
face and the nodes of the gap intersect the Fermi surfaces 
near x and y axes. The peak position of the spin suscep- 
tibility Xs{q) shown in Fig. 4, which should correspond 
to the nesting vector of the Fermi surface, stays around 
Q = (w, 7r/4) regardless of the pressure. This vector Q 
bridges the portion of the Fermi surface with <p(k) < 
and (j)(k + Q) > 0, which is the origin of the d xy -~w&ve 
like gap. 
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(a) 14.1GPa 5.8KU=0.7eV 



(b)16.8GPa7.2KU=1.2eV 
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Fig. 4. Contour plots of the spin susceptibility Xs(<j) f° r the same 
values of U and p as in Fig. 3. 
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Fig. 5. Contour plots of |G(fc)| 2 , the nodes of the gap function, 
and the spin susceptibility at (a,)U = 0.7 eV, p = 14.1 GPa. 
(b)U = 1.2 cV,p = 17.2 GPa. At the circled portions of the 
Fermi surface, (f>(k)ip{k + Q) is positive. 



4. Discussion 

In this section, we physically interpret some of our 
calculation results. 

4-.1 Pressure dependence of T c at fixed values of U 

For fixed values of U, T c tends to be suppressed at high 
pressure as seen in Fig. 2. This can be explained as fol- 
lows. We have seen in Fig. 4 that the peak position of the 
spin susceptibility does not depend on pressure, but the 
peak value itself decreases with pressure for fixed values 
of U as shown in Fig. 6. This is because the nesting of 
the Fermi surface becomes degraded due to the dimen- 
sional crossover of the Fermi surface mentioned previ- 
ously. Consequently, the pairing interaction (nearly 
proportional to the spin susceptibility) becomes smaller, 
so that T c becomes lower, with increasing pressure. 

For U = 0.9 eV (and possibly for U = 0.8 eV), T c is 
slightly suppressed also at low pressure, so that a opti- 
mum pressure exists. This may be due to the fact that 
at low pressure, Q (spin susceptibility peak position) 
bridges some portions of the Fermi surface that has the 
same sign of the gap (Fig. 3(a)). In fact, this tendency of 
Q bridging the same gap sign is found to be even stronger 
for lower values of pressure as seen in Fig. 5(a), while at 
higher pressure, the nodes of the gap run along the Fermi 
surface so as to suppress this tendency (Fig. 5(b)). We 
will come back to this point in §4.3. 

4-2 Pressure dependence of the maximum T c upon vary- 
ing U 

As can be seen from Fig. 2, for each value of pressure, 
there exists an optimum value of U(= U op t) which max- 
imizes T c . In this subsection, let us discuss the pressure 
dependence of this optimized T c as a function of Z7 opt , 
that is, T c ([/ opt ). Regarding the lower pressure region, 
Xo is large because of the good nesting of the Fermi sur- 



face, so that U has to be small in order to avoid AF 
ordering (in the sense mentioned in §2). This is the rea- 
son why the maximum value of T c is relatively low in 
the low pressure regime as seen in Fig. 2. On the other 
hand, in the high pressure region, xo is small because of 
the 2D-like Fermi surfaces, so that U must be large in 
order to have large \ an d thus pairing interaction. Such a 
large U, however, makes the normal self energy S(fe,£„) 
large, which again results in low T c (In Fig. 5.(b),the low 
height of \G(k)\ 2 represents the large effective mass of 
Fermion.). Thus, relatively high T c (U opt ) is obtained at 
some intermediate values of pressure. 

4-3 Pressure dependence of the gap function 

In this subsection we discuss the variation of the gap 
function with increasing pressure. To understand this 
variation in real space, we use the following relation 

O = ^0(fc)c feT c_ fei 
k 

= ^29(8){CilCi+5l ~ CiXQ+sj), (11) 
i,S 

where i and i+5 denote sites in real space where pairs are 
formed, and g(S) is a weight of such pairing. Note that a 
'site' here corresponds to a unit cell (or a dimer). Con- 
sidering up to 22 nd nearest neighbor pairings, we have 
determined a set of g(S) that well reproduces <j)(k) ob- 
tained by the FLEX, using least squares fitting, as shown 
typically in Fig. 7. 

The result of this analysis is shown in Fig. 7, where the 
thickness of the lines represents the weight of the pair- 
ing determined from the value of g(S). We can see that 
the direction in which the dominant pairings take place 
changes from b to b+c as the pressure increases, which 
looks more like d xy -wa,ve like pairing. These changes of 
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Fig. 6. Pressure dependence of the 
Xs(q,u) m = 0) for several values of U. 
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the pairing directions are correlated with the increas- 
ing of hopping t(p2). Thus, from the viewpoint of this 
real space analysis, we can say that the change of the 
dominant pairing direction due to the increase of t{p2) 
suppress the tendency of the nesting vector Q bridging 
the portions of Fermi surface with the same gap sign. 

4-4 Origin of the "high T c " 

Finally, we discuss the reason why the obtained results, 
namely the values of T c and the form of the gap function, 
resemble that of the dimer limit approach. The present 
situation is in sharp contrast with the case of k-(BEDT- 
TTF)2X, where it has been known that, compared to the 
results of the dimer limit approach, 14-16 the position of 
the gap nodes changes and the values of T c , if any, is 
drastically reduced in the original four band model with 
moderate dimerization. 8 ' 9 

A large difference between the present case and k- 
(BEDT-TTF) 2 X is the Fermi surface nesting. As men- 
tioned previously, the quasi-one-dimensionality of the 
system gives good Fermi surface nesting with strong 
spin fluctuations, fixing the nesting vector and thus the 
pairing symmetry firmly, while in the case of k-(BEDT- 
TTF)2X, the Fermi surface has no good Fermi surface 
nesting. 

However, it is unlikely that the good Fermi surface 
nesting alone can account for the resemblance between 
the 3/4- filled model and the dimer limit model because, 
for example, in the study for another organic super- 
conductor (TMTSF) 2 X, it has been known that a 1/4- 
filled model with no dimerization gives weak spin fluc- 
tuations within FLEX even though the Fermi surface 
nesting is very good. 17 One difference from the case of 
(TMTSF^-X' is the presence of moderate (not so weak) 
dimerization, but there is also a peculiar structure in the 
density of states as pointed out in the previous study. 4 
Fig. 8(a) shows the density of states of the antibonding 
band at p = 17 GPa for [7 = 0. The two peaks near the 
Fermi level (energy=0) originates from the saddle points 
of the band dispersion located at the T point and the 
Y point (k = (-tt, 0)). Consequently, the "Fermi surface 
with finite thickness" , defined by Ep — S E < E(k x ,k y ) < 
Ep + SE, becomes thick near the T and the Y points, as 
shown in Fig. 8(b). In fact, this trend is already seen in 
the contour plots of the Green's function in Figs. 3 and 5, 



Fig. 7. (upper panel) Results of the least squares fitting of the gap 
function. The dashed curve is the result of the FLEX calculation 
and the solid curve is the fitting, (lower panel) The weight of the 
directions in which the pairings take place for (a)p = 14.2 GPa, 
U = 0.7 eV (b) p = 17.2 GPa, U = 1.2 eV 



where the bundles of the contour lines become thick near 
the r and/or the Y points. Importantly, the wave vector 
(the nesting vector ~ (tt, tt/4)) at which the spin fluctua- 
tions strongly develop bridges the states near the T point 
and those somewhat close to the Y point (Fig. 8(b)), so 
that there are many states which contribute to the pair 
scattering. From the above argument, our results suggest 
that the coexistence of the good Fermi surface nesting, 
the large density of states near the Fermi level, and the 
moderate dimerization cooperatively enhances electron 
correlation effects, thereby giving results similar to those 
in the dimer (strong correlation) limit. 

Now, these factors that enhance electron correlation 
should also make T c itself rather high. In fact, T c of 
~ 0.0006W^ almost reached in the present model, where 
W is the band width (around 1.3 eV for p = 16 GPa), is 
relatively high among T c obtained by FLEX+Eliashberg 
equation approach in various Hubbard- type models. 
Namely, Arita et al. have previously shown 18 that T c of 
order O.OOIH^ is about the highest we can reach within 
the Hubbard model, 19 which is realized on a two dimen- 
sional square lattice near half filling, namely, a model for 
the high T c cuprates. The present study is in fact remi- 
niscent of the FLEX study of the high T c cuprates, where 
the 3/4- filled two band model 20 and the half-filled single 
band model indeed give similar results on the supercon- 
ducting T c and the pairing symmetry. 3 The cuprates also 
have a large density of states at the Fermi level originat- 
ing from the so called hot spots around (tt, 0) and (0, 7r), 
and the wave vector ~ (ir, it) at which the spin fluctua- 
tions develop bridges these hot spots, as shown in Fig. 
8(c). Moreover, a moderate band gap also exists in the 
cuprates between the fully filled bonding/non-bonding 
bands and the nearly half-filled antibonding band. The 
situation is thus somewhat similar to the present case. 
To conclude this section, it is highly likely that the co- 
existence of the factors that enhance correlation effects 
and thus make the results between the 3/4- filled original 
model and the half-filled dimer model similar is the very 
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Fig. 8. (a) The density of states of the antibonding band at p = 
17 GPa for (7 = 0. (b) The Fermi surface with finite thickness 
defined by Ep — SE < E{k x ,k y ) < Ep + 5Ep, where E(k x ,k y ) 
is the band dispersion, and SEp = 0.015 eV is taken here, (c) A 
similar plot for the cuprates. 



reason for the "high T c " in the title material. 

5. Conclusion 

In the present paper, we have studied the pressure 
dependence of the superconducting transition tempera- 
ture of an organic superconductor /3'-(BEDT-TTF) 2 ICl 2 
by applying two-band version of FLEX to the original 
two-band Hubbard model at 3/4-filling with the hop- 
ping parameters determined from first principles calcu- 
lation. The good Fermi surface nesting, the large density 
of states, and the moderate dimerization cooperatively 
enhance electron correlation effects, thereby leading to 
results similar to those in the dimer limit. We conclude 
that these factors that enhance electron correlation is the 
origin of the high T c in the title material. 

As for the discrepancy between the present result and 
the experiment concerning the pressure regime where the 
superconducting phase appears, one reason may be due 
to the fact that we obtain T c only when Uxo > 0.995 is 
not satisfied despite the fact that this criterion for "anti- 
ferromagnetic ordering" , originally adopted in the dimer 
model approach, 4 does not have a strict quantitative ba- 
sis. Therefore, it may be possible to adopt, for example, 
Uxo > 0.999 as a criterion for "antiferromagnetic order- 
ing" , which will extend the superconducting phase into 
the lower pressure regime. Nevertheless, it seems that 



such consideration alone cannot account for the discrep- 
ancy because in order to "wipe out" the superconduc- 
tivity in the high pressure regime as in the experiments, 
smaller values of U would be necessary, which would give 
unrealistically low values of T c . Another possibility for 
the origin of the discrepancy may be due to the ambigu- 
ity in determining the hopping integrals from first prin- 
ciples calculation. Further qualitative discussion may be 
necessary on this point in the future study. 
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